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1. Introduction 

Let (M, g, 51) be the conformal rescaling of a 3+1 dimensional vacuum spacetime 
(M,g), with g = ri^g, where 51 is a smooth function on M vanishing on the 
boundary 9M. If M has smooth null boundary J^, with connected components of 
topology S'^ X M, it follows that the Weyl tensor of (M, g) vanishes on . From this 
follows peeling properties for (M,g). This picture of isolated systems in general 
relativity, developed by Penrose, see jT6[ |9|, is useful for studying the mass and 
angular momentum of spacetimes, as well as gravitational radiation. 

Friedrich has found a first order symmetric-hyperbolic version of the Einstein 
evolution equations, called the conformally regular field equations, which may be 
extended through ^ ||l^. The Cauchy data for this system on a future Cauchy 
surface M, asymptotic to so that M n 9M = dM, include gab, Kab, 51 as well 
as components of the rescaled Weyl tensor 51~^C"^^^. In order to get a regular 
evolution at J^, these data must be regular up to dM. It was shown by the author, 
Chrusciel and Friedrich H, ^, |^ that under certain conditions on the boundary 
geometry of M, the Cauchy data for the conformally regular field equations, are 
smooth up to dM. 

Using the conformally regular field equations, Friedrich has shown that the max- 
imal vacuum development of data (called hyperboloidal data) on a future Cauchy 
surface intersecting J^, regular up to dM is a spacetime which has a "smooth 
piece of J^" . Further, for small data the maximal vacuum development has a null 
boundary with future complete null generators and a regular timelike infinity. 

A programme has been initiated to numerically evolve the Einstein equations us- 
ing the conformally regular field equations ^ |l3, 12 . This approach may have 



advantages over the "traditional" approach which treats the asymptotic flatness 
condition by introducing boundary conditions far away from the isolated system 
under study, and attempts to observe for example gravitational wave signatures on 
this boundary. 

In this note we will discuss the conformal procedure for constructing solutions 
{gab, Kab, 51) to the constraint equations and the geometric conditions for regularity 
at J^. I will also briefly discuss the Cauchy problem for the Einstein evolution 
equations at J^, directly from the point of view of the Einstein equations in (M, g). 



This paper is based on a talk given at the Workshop on the Conformal Structure of Space- 
Times held at Tiibingen, April 2—4, 2001. I am grateful to the organizers, Jorg Prauendiener and 
Helmut Friedrich, for their hospitality and support. This work is supported in part by the Swedish 
Natural Sciences Research Council (SNSRC), contract no. R-RA 4873-307 and NSF, contract no. 
DMS 0104402. 
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2. Preliminaries 

The following index conventions will be used. Greek indices a, (3, . . . run over 
0,...,3, lower case latin indices a,6, c, ... run over 1,2,3, and upper case latin 
indices A,B,C,... run over 2, 3. 

Let (M,g) be a 3+1 dimensional globally hyperbolic space-time with covariant 
derivative D. Introduce coordinates (x") = {t, x") on M, and consider the foliation 
Mt consisting of level sets of t. We will often drop the subscript t on Mt and 
associated fields. M will be a compact manifold with boundary dM, and closure 
M = M U dM. Let T be the normal of M, assumed to be timelike, and let g be 
the induced positive definite metric on M, with covariant derivative V. We will use 
an index T to denote contraction with T, for example It = T'^ta. Let S = D^ri 
and assume jSj > 0. The lapse function N and shift vectorfield X are defined by 
dt — NT + X. The second fundamental form K is defined by Kab — (DaebjT). 
Then we have Kab = -^Crgab = -^N'^^dt - X)gab and DaCfc = Vaef, - TKab- 

Assume that M has null boundary dM and let M = M U dM. Let n he a 
positive function on M with SI = 0, d^l ^ on dM and denote by g the conformally 
related metric g = 0~^g. We will call g the unphysical metric and g the physical 
metric. We will use (•,•) for the inner product induced on TM by g. Geometric 
quantities associated with g will be decorated with a tilde ~, for example the 
covariant derivative D and the Ricci tensor Rq,^. We will consider only the case 
when (M, g) satisfies the vacuum Einstein equations Rafi = 0. Let uj — 
Assume that dM = M O dM, has positive definite induced metric and that ui is 
a defining function for dM, i.e. duj\gj^,^ ^ 0. We will refer to dM as and the 
surface dM C dM as a cross section of J^. 

There is a gauge ambiguity in the conformal rescaling. Let Ct — i76~^, g — O^^g 
for some positive function which is bounded and bounded away from zero on M. 
Then g — Ct^^g is another conformal rescaling of g. A conformal gauge change 
can be used to control the mean curvature of dM , as well as the mean curvature 
of M. In particular, the conformal gauge freedom can be used so that the level 
sets of uj are the leaves of the gauss foliation w.r.t. dM, i.e. there are coordinates 
(cc") = {x,y^), so that in a neighborhood of dM, we have ui = x, and the metric 
takes the form 

(1) gabdx°-dx^ = dx^ + hAsdy^dy^ . 

This possible without loss of generality, cf. H, Lemma 2.1]. By construction we 
have x{p) = d{p, dM) and the unit normal dM is rj = dx- We use the index 1 for 
contraction with rj, for example tu ~ \7°-xtab — v'^tab- Let y a and fiAB denote 
the covariant derivative with respect to Hab and the Ricci tensor of hAs- The 
second fundamental form of dM with respect to rj is defined by A as = (Vyies, ry). 
Then we have V^es = y a^b + V^ab and Xab = —^^r^hAB- 

The quantities associated to {gij , K^j , N, T) in the physical space-time are 

(2a) (jij = uj'^gij, Kij = uj^^Kij + ^gij, 

(2b) N = Lj-^N, f = ujT. 

The shift vectorfield X does not scale, so X — X . The physical Weyl tensor satisfies 
Cap-iS = R.a/375 siucc Ra/3 = 0. Further, C"^^^ = ^""p-yS- 



CONSTRUCTION OF HYPERBOLOIDAL INITIAL DATA 



3 



A symmetric tensor field Aat on M is said to satisfy the shear free condition if 



(3) {Aab ~ i/i^^AcD/iAB)|aM = ^M^b - ^^''XcDhAB)] 



dM- 



In the foUowing we will consider only hypersurfaces intersecting J^_|_ which in view 
of our conventions force S < 0, so that — —1. The spacetime (M, g) is said 

to be shear free at dM if Kab satisfies the shear free condition. In order for g to 
be in C^(M), it is necessary that Kab satisfy the shear free condition, cf. Prop. 
3.1]. In the following we only consider initial data which are shear free. Note that 
a change of conformal gauge can be used to get h^^^ \ab\qj^[ — 0. 

3. Conformal rescalings of Minkowski space 

In order to understand the Einstein equations near ^ it is important to choose 
an appropriate conformal compactification and a suitable foliation near ^ . In this 
section we will display a few examples in the simplest case, Minkowski space. 

Minkowski space is M'' with the flat line element 

ds^ ^ -dt^ + dr"^ + r^da^, 

in radial coordinates {t,r,9,(/)), where da^ — d6^ + sin^ {9)dcj)'^ is the line element 
on 5^. We will consider two conformal rescalings of Minkowski space. First, define 
new coordinates r, p by 



t + r — tan — - — , t — r — tan 



T - p 



(4) U ^ {t,p : -TT < T < TT, < p < TT, \t + p\ < TT , |r-p|<7r}. 
Rescaling with the conformal factor 

J7(r,p) = 2cos(^^)cos(^^ 

which is positive on U, gives 0, p. 450-451] 

ds^ = gai3dx°'dx'^ = -dr^ + dp'^ + sin^(p)da-^, 

which is the line element of the Einstein universe M x 5'^. We have N = 1, X = 0. 

The null part of the boundary of U has two components with t > or r < 0, 
called future and past null infinity, denoted J^±. The point (t, p) = (0,7r) is called 
spatial infinity and denoted zq. The points i± with (r, p) = (±7r, 0) are called future 
and past timclike infinity. All inextendible future directed null curves start at 
and end at ^+ . 

The hyperboloids — = in Minkowski space correspond to hypersurfaces 
in U intersecting dU at t — tt/2 at different angles depending on u, cf. figure ^ 
In particular, the hyperboloid with sectional curvature —1 corresponds to r = tt/2. 
On the other hand, the time translated hyperboloids given hy t = vl + r^ + u 
intersect at r = ■k/2 + arctan(u), see figure |[ In this coordinate system, the 
future translated hyperboloids approach timelike infinity i+ as u increases. The 
Lapse N is given at by 

lim N = ^ = cos(r - 7r/2) 



which tends to zero as r increases to 



TT. 
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Figure 1 . Hyperboloids — r^ = v? in conformally compactified 
Minkowski space. 
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Figure 2. Time translated hyperboloids t = V 1 + + u in con- 
formally compactified Minkowski space. 

For the second conformal rescaling, introduce new coordinates if}, r by 

t = coth{%li) + T, 
r = 1/ sinh('(/)). 

This is a singular coordinate transformation, but this does not cause a problem as 
we are only interested in the behavior near infinity. In terms of the new coordinates, 
the Minkowski line element takes the form 

df- = -ch^ + 2 sinh-2(V))dVrfT + sinh-=^(V')(#^ + dcr^). 

After a conformal rescaling with 

Q, = smh(ip), 

we get the line element ds^ — fPds^ given by 

(5) ds^ = ~ sinh^ {^)dT^ + 2d^dT + (dV/ +da'^). 

The hypersurfaces r =constant correspond to the time translated hyperboloids 
t = Vl + r"^ + T in Minkowski space. The conformally rescaled spacetime M has 
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Figure 3. Hyperboloids 



in the conformal rescahng (|5|). 



null boundary ^+ — {ip = 0}, and so the boundary of this conformal rescaling does 
not include or io, i±- Further, the closure of M is not compact. 

The foliation AIj- of level sets of r has Lapse N — cosh('(/;), and Shift X — 
so that the time like normal T is T = cosh~^(V')(i9T — 9^)- The foliation Mr is 
static so that Kab = 0. The induced geometry on Mr is the cylinder c?-0^ + dcr^, 
with totally geodesic boundary = 0}. The future sheets of the hyperboloids 
<2 _ 7.2 _ y2 Q-j- j-Qga^j]^ curvature — 3/u intersect at t = 0, see figure |[ 

4. Conformal constraint equations 
The Einstein vacuum equations Rq/3 = imply the constraint equations 
(6a) R+{k\f -kabk'^'- ^Q, 

(6b) V"i^ah - ^bk'c = 0- 

We will consider the constraint equations, in the special case of constant mean 
curvature Votri^ = 0, as a system of equations for the unphysical data and rewrite 
them in a form which exhibits the freely specifiable data. 

In the form of the constraint equations that we will present, the freely specifiable 
data are ((7^,^^), where gij is a metric on M, and A^j is a trace free, symmetric 
tensor field on M satisfying the shear free condition. We will assume for simplicity 



dM 



0, which does not restrict the degrees of freedom. 



Let p be a defining function for DM, p = x near dM. Assume Ej^^^ is a nonzero 

constant. Set gab — p^^gab, define V, R etc. w.r.t. g and raise and lower indices on 
hatted fields with g. Given data gij,Aij, a solution {gij, Kij,uj) to the constraint 
equations is constructed by solving the system (|^) below. We will refer to the system 
(^ as the conformal constraint equations. Existence, uniqueness and regularity of 
solutions to the conformal constraint equations has been analyzed in the papers 
|j. For a vector field Y, let 



(7a) 



For a solution Y to the system 

(7b) V^L(F),, 
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define aij by 

(7c) cTy = - p^^L{Y),j. 

Then aij — p~^aij satisfies V"(Tah — 0. Let u be tlie solution to tlie semilinear 
elliptic equation 

(7d) -8Au + iku- a,ja'^u-'^ + GS^u^ = 0. 

Finally, define {gij, Kij,uj) hy lo = u^'^p, Kij — u~'^aij + 'Sgij. The fields {gij, Kij) 
corresponding to a solution gij , Kij , oj of the conformal constraint equations (|7|) 
solve the physical constraint equations (|^) with mean curvature g^-' Kij = 3S, under 
the constant mean curvature condition VqS = 0. 

4.1. Constant mean curvature hypersurfaces. The existence of complete con- 
stant mean curvature hypersurfaces in asymptotically flat spacetimes was studied 
in [0. There it was proved that in an asymptotically Schwarzschild spacetime 
satisfying a natural causal condition (domain of dependece of compact sets is com- 
pact), given a cross section of and a number E < 0, there is a unique CMC 
hypersurface intersecting ^+ at the given cross section. The proof uses a barrier 
construction, together with the causality condition, to prove that the solutions to a 
sequence of boundary value problems converges to a complete CMC hypersurface. 

This result indicates that a wide class of spacetimes with regular conformal 
compactification are foliated by CMC hypersurfaces. The regularity of the height 
functions of these CMC hypersurfaces at ^ has not been studied in detail, nor has 
the existence proof been carried out in the conformally compactified setting. This 
is a natural problem which should be studied further. 

The conditions for regularity of CMC solutions of the conformal constraint equa- 
tions constructed using the procedure discussed above is well understood, see the 



discussion in section 4.3 below. However, it is not at all well understood how these 
results relate to the regularity of CMC foliations in the physical or unphysical 
spacetime, evolved from these data, nor what spatial gauge conditions are suitable 
in this situation, see section |^ below. 

4.2. Degenerate elliptic equations. The system (^ of conformal constraint 
equations is a degenerate elliptic system for (u,y). We will briefly describe the 
regularity properties of solutions to systems of this type. As an example consider 
the Lichnerowicz equation ([zd]). The principal part of the right hand side of ( [z^ ) 
is essentially of the form 

L = L^ + B, 

where is the ordinary differential operator = x'^d'^ + axdx + b, and B is of 
lower order in dx- In case (1 — a)^ > 45, the equation for L^x" = has real distinct 
characteristic roots a± and the equation L^u — f has a solution u = o{x°'^) for 
sufRciently regular / — o{x°'^ ). 

The operator appearing in the Lichnerowicz equation has critical exponents a- = 
— 1, a+ = 3. Let L be this operator and consider the equation Lu — /, for 
/ e C°°(M). Suppose / is of the form / = fix + f2x'^ + fj,x^. Then a solution 
u = o{x^^) is of the form 

u — 1 + uix + U2X^ + u^,ix'^ Inx -|- u^x^ + higher order terms, 

where W3.1 = c/3 for some explicit constant c. This example reflects the fact, 
cf. H, |[ ^ that solutions to degenerate elliptic systems are in general nonsmooth 
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at dM , instead the general form of the solution has a poly logarithmic expansion of 
the form u — 'Y^ u^jX^ In-* x. In case the RHS and the coefficients are smooth, the 
logarithm terms in u appear first at the critical exponent a+. 

The terms in the expansion of u up to and including the first logarithm term are 
computable in terms of the data, which implies that explicit geometric conditions 
for the regularity of solutions of the conformal constraint equations , and hence 
of the physical spacetime constructed from these data, can be found. 

4.3. Regularity of solutions to the conformal constraint equations. We will 
now discuss the detailed conditions for regularity of the conformal data {gab^ Kab, ^) 
at dM. We will assume that free data {gij,Aij) G C°°{M) are given and consider 
the conditions needed for the solution (fc, Kij,uj) of the conformal constraint equa- 



'^AB — ^AB — t;^^^ ^CD^AB 



tions, produced by solving the system ([7|), to be in C°°{M). 

In the case when M is a moment of time symmetry in the unphysical space-time, 
i.e. Aij = Kij — 0, there is a simple criterion for regularity of the solution to the 
constraint equations. Let 

1 

2' 

denote the trace free part of \ab- Then the solution to the conformal constraint 
equation is in C°°{M) if and only if the conformal density C on 9Af , defined by 

c = y ^Xab + A*-4Si?^5 _ \h^''\cD\^''XAB 

vanishes. In particular, this holds if DM is totally umbilic, i.e. if 

(Aab - -j^h'^'^ ^cd^ab)\qj^j = 0. 

In general if this condition does not hold, then the solution has an expansion in 
powers of x and \nx near dM and is thus of finite regularity at dM. 

In the case where Kij ^ and consequently Aij ^ 0, we make the following 
simplifying assumptions. We consider shear free data {gij^Aij), such that 

(8a) /i^^AAB|g^,, = 0, 

(8b) ^i.leM = 0- 

Equation should be viewed as a conformal gauge condition, which can always 
be satisfied if the data (g^ , Aij ) is smooth up to boundary. Similarly, equation 
( pb| ) may be viewed as a gauge type condition as it is possible to achieve this by 
subtracting a term p^^L^ from Aij. If (|8|) holds, then the fields uj^Y^ satisfy 

(9a) Lo ^\Y.\x + 0{x^), 

(9b) = 0(a;3), 

and Kab is determined up to first order by Aab- If (||) holds, then the fields 
oj,Kij,LO~^Cai3'yS are in C°°{M) if and only if the equations 

(10a) O^Kab - \h''''d^KcDhAB = 0, 

(10b) y-^y^A^B+i^AsA^^-O, 

(10c) yid^KAB + \\cD\''''hAB) = 0, 
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cf. |, Eq. (4.21)-(4.23), (4.33)], hold on dM. Note that our sign conventions for 
e and Kab are opposite those of M. The particular form we present here holds 



assuming (ph. Equation (10a), which states that OxKab is proportional to Hab on 



9M, is a consequence of the vanishing of the Weyl tensor on dM, while equations 



(10b) and (tOd) follow from smoothness of Kab up to dM. In particular, it can be 



shown that if dM = S^, then 
(11) (^d^KAB + ^A^^Acu/iab) la^, = 0. 

5. The initial value problem 

In this section we will consider the Einstein evolution equations directly from 
the point of view of the unphysical Cauchy data. Since we know from the work of 
Friedrich [0 that the maximal vacuum extension of conformally regular data on M 
has regular conformal boundary near dM, it follows that solutions to a well-posed 
formulation of the Einstein evolution equations in the unphysical space-time will 
have this property. Working directly in terms of the unphysical Cauchy data allows 
one to relate this problem to the extensive literature on the numerical solution of 
the Einstein evolution equations. In particular, it is interesting to consider gauge 
choices and hyperbolic reformulations. 

Consider a foliation Mt of (M, g) with dt = NT + X. The structure equations 
of the foliation imply the Einstein evolution equations 

(12a) dtg^J = ~2NK,j + CxQrj , 

(12b) dtK^j = -V,Vj7V + N[R,j + tvKK^j - 2K„nK"; - Ry) + CxK^J, 
and the constraint equations 

(13a) R~\K\^ + {tTKf ^2IiTT + 'R, 

(13b) V^trK - V^Xy = Rn. 

The unphysical Ricci tensor R^^ is given by 



R„^ = R^^ - n-^[2T>a'Dpn + 'D-^'D'^nsafs] + 'sn^^n^nD'^nga 



/3> 



the unphysical Ricci scalar is R = g^'^Ra/j, and — KijK^^ . The condition 
T)t^ = S may be viewed as an evolution equation for f2, 

(13c) dfO. = NT. + XVl. 



5.1. Gauge condition at dM. It is natural to require 9til|g^^ = 0. From (13c) 
and uj = \T\x + h.o., 9t^i|^J^^^ = iVS + so 9tri|g^_j = impHes 

In view of this, a natural boundary condition for N and X is 
(14) = ^laM = ^' 

which corresponds to the asymptotic behavior of Lapse and Shift in the second 
conformal compactification discussed in section ^. In a neighborhood of dM, we 
may decompose X as X = arj + (3. Then the boundary condition for X is a = 
l,/3 = on dM. 
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5.2. Evolution at 9M. As discussed in §5] , the shear free condition is necessary 
in order for the development of the data on M to have a regular conformal boundary. 
It is convenient to introduce the notation A = h'~'^XcD, k = h'^^KcD- We relax 



the condition A|g^^ = used in section |4. 



Imposing the boundary gauge condition (|14[), a calculation shows that 

dthAB\gj^j = + X)hAB\gj^j- 

In particular, dthAB\gj^f is pure trace. For a symmetric tensor tAB, let — 
tAB — ^h^'^tcohAB denote the trace free part. Define the shear tensor aAB by 

a-AB = X*j^B + ^AB, 
so that the shear free condition can be formulated as 

From the fact that dthAB\gj^j is pure trace together with the fact that M is 3- 
dimensional, it follows that 

^t<JAB\gJ^,J = [dt{KAB + Xab)]*- 

Thus, if dt{KAB + -^ab) Ig^^ is pure trace, then the shear free condition is conserved 
by the evolution. The evolution equations and the boundary conditions imply 

dt\AB\gj^,j = dxi^AB + Kab) + NiKab - ^Cxi9AB\gj^f, 

where A^i = (9x-^|gj,j, Xi = X] j^^^^. After a lengthy calculation one finds, 

dtCTABlgj^j = NiK\b - ]^[Cx^9ab]*- 

The leading order term Ni is determined by the regularity of Tt^p up to (9M. 
Therefore the condition for preserving the shear free condition under evolution, 
StCTyis Ig^j ~ shows up as a gauge condition for X. Due to the shear free condition, 
the equation 

NiK\b - 2^^x,9AB]*\gj^j = 0, 

can be solved by putting Xij^^^ — Nii], in other words by imposing the condition 
that X is parallel to 77 at dM to first order, and that dt is null to first order. 

6. Discussion 

In this note, I have shown how to construct solutions of the conformal con- 
straint equations, and indicated the first steps in analyzing the evolution problem 
for the Einstein vacuum evolution equations, directly in the unphysical setting. It 
seems worthwhile to explore gauge choices, hyperbolic reformulations, boundary 
conditions etc. for the evolution equations directly in this setting. This makes it 
possible to tie in to the extensive literature on the standard form of the Einstein 
evolution equations. As an example, it is interesting to look at reformulations of the 
Einstein equations, derived by adding multiples of the constraints in the evolution 
equations, which have improved stability of the constraints, see for example It 
is natural to ask here if there are such reformulations which are also better behaved 
at J^. 

The gauge conditions that have been considered so far in the literature on the 
numerical solution of the Einstein equations are mainly hyperbolic in nature, and 
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therefore suffer from well known problems such as gauge singularities It 
is known that CMC time gauge gives a well posed form of the Einstein evolution 
equations. Therefore it seems natural to consider this time gauge also for the 
hyperboloidal IVP. 

Choosing a time gauge condition such as the CMC condition, leaves open the 
choice of spatial gauge. It seems likely that it is advantageous to work with a "fully 
gauge fixed formulation" for numerical work. The combination of CMC time gauge 
and harmonic spatial gauge gives a well posed elliptic-hyperbolic system in case of 
compact Cauchy surface Does this gauge condition give a well posed system 
also for the hyperboloidal IVP, with boundary conditions given by a foliation of 
J'l The defining equation for the Shift vector field given in is derived from 
the evolution equation for the metric and therefore does not depend on the matter 
content of the spacetime. This makes it possible to apply this formulation both in 
the physical and unphysical spacetime. It is an interesting problem to explore the 
consequences of this gauge choice for the hyperboloidal IVP, both from the point of 
view of the physical and unphysical formulation of the IVP and gauge conditions. 
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